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<^ • Abstract 

We consider a <i-dimensional random quantum walk with site-dependent ran- 
dom coin operators. The corresponding transition coefficients are characterized by 
1 deterministic amplitudes times independent identically distributed site-dependent 

random phases. When the deterministic transition amplitudes are close enough 
to those of a quantum walk which forbids propagation, we prove that dynamical 
localization holds for almost all random phases. This instance of Anderson local- 
ization implies that all quantum mechanical moments of the position operator are 
uniformly bounded in time and that spectral localization holds, almost surely. 
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Quantum walks denned on a lattice, deterministic or random, have become a popular 
topic of research recently, due the interest they have for different scientific communities, 
see the reviews [26], [29], [3]. Quantum walks can be used to provide an effective de- 
scription of the dynamics of certain types of physical quantum systems and belong, in 
this sense, to the broader category of quantum network models, [1], [33], [13], [9]. In 
particular, quantum walks accurately describe the dynamics of atoms trapped in certain 
time dependent optical lattices, as was recently demonstrated experimentally, [24], [41]. 
Quantum walks also offer a field of investigation for probabilists since the quantum in- 
terpretation of the wave function together with the relative simplicity of the dynamics of 
quantum walkers give rise to generalizations of classical random walks possessing some- 
times unfamiliar features, like ballistic transport properties, see e.g. [29], [8]. Moreover, 
in the same way as classical random walks are important for classical computer science, 
quantum walks play an important role in the benchmarking procedures used to assess 
quantum algorithms in theoretical quantum computing [35], [32]. 

It is of interest in some of the instances mentioned above to introduce static random- 
ness in the dynamics of the quantum walk, in order to take into account fluctuations 
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in the model, be they unavoidable experimentally or inherent to the model, [34], [40]. 
In such cases, one expects the transport of the system to be completely suppressed due 
to the destructive interferences induced by the randomness and one speaks of localiza- 
tion. This well known phenomenon in solid state physics goes under the name Anderson 
localization and it is mathematically well established for a large variety of continuous 
time dynamics driven by random Schrodinger operators, in various regimes, see e.g. [39], 
[27]. While of a similar nature, this phenomenon has been less studied for the case of 
discrete quantum dynamics defined by random unitary operators. Localization results 
have been proven for certain types of random unitary operator in [10] [19], [38], [16], [17], 
[6] , [7] with various physical and mathematical motivations, see [20] for a review. Con- 
cerning coined random quantum walks specifically, the following contributions address 
such questions, [28], [22], [37], [2]. The first paper yields a special instance where the 
randomness fails to induce localization, while the last three papers provide localization 
results. These results hold for the dimension one only. 

Our purpose in this paper is to prove dynamical localization for d- dimensional ran- 
dom quantum walks, d > 1, characterized by random transition coefficients and de- 
terministic transition probabilities, in the same spirit as the one-dimensional model 
addressed in [22]. We make use of ground work on dynamical localization by means of 
fractional moments estimates for random unitary operators provided in [17] and of the 
methods used in the recent paper [7] to tackle two dimensional problems in a perturba- 
tive regime analogous to the large disorder regime of the Anderson model. 

Let us close this introduction by briefly mentioning two unrelated instances where 
one speaks of localization and of randomization in the study of quantum walks and 
which are not discussed in this paper. First, under some circumstances, deterministic 
quantum walks that are either homogeneous in space or spatially locally perturbed by 
the presence of defects give rise to dynamics with the following features: An initial states 
is split into a part that travels to infinity ballistically and a part that remains for all 
times in a neighborhood of the origin. The persistence of a part of the quantum walker 
close to the origin, sometimes also called localization, is of a completely different nature 
from Anderson localization and is due to the appearance of eigenvalues in the spectrum 
of a relevant unitary operator, in addition to the usual continuous spectrum, see [18], 
[10], [11], [30]. The second point which is not discussed here concerns a different way of 
randomizing a coined quantum walk. It consists in making the coin matrices depend on 
time in a random way, which leads to the study of random dynamical systems for the 
corresponding random quantum walk. Various models of this sort have been considered 
in the physics literature and numerical as well as rigorous studies have been lead, see 
e.g. [36], [31], [3], [21], [15], and references therein. The main outcome is the appearance 
of diffusion phenomena in the averaged dynamics, restoring a flavor of classical walk. 
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2 The Model and Main Result 



We consider random quantum walks on the regular lattice 7L d with coin space C 2d , i.e. a 
unitary random operator on % = C 2d (g) l 2 (Z d ) ~ l 2 (Z d ; C 2d ), by means of the following 
construction. 

Let us denote the standard basis of C 2d by 

{|r)} re/± , where I± = {+1, -1, +2, -2, . . . , d, -d} (2.1) 

and that of l 2 (7L d ) by {\y)} ye z d - The corresponding basis vectors of the tensor product % 
will be written as \r,y) = \r) ® \y) and orthogonal projectors on these basis vectors are 
denoted by P T = \t)(t\ etc... The scalar product of tp, ijj G % is denoted by (ip,xjj), with 
linearity on the right hand side, and the associated norm is denoted by ||y?|| 2 = (<f,<f)- 
We also denote by || ■ || the norm of bounded operator on H. We set r : I± — > Z d by 

r(r) =sign(r)e w , r G 7±, (2.2) 

with {ej}j =1: 2,... : d the canonical basis of R d . 

Let S* be the coin state dependent shift operator defined on V, by 

S= Yl Pr®\y + r{r))(y\. (2.3) 

yeZ d ,r£l± 

In other words, S\a,x) = \a, x + r(a)), for any basis vector \a, x). Let C G U(2d), U(2d) 
the set of unitary matrices on C 2d , be a coin matrix. Then the unitary operator § on "H 
given by 

§ = S(C®I)= Yl PrC(S)\y + r(r))(y\ (2.4) 

yei, d ,rei± 

defines a deterministic quantum walk. 

Consider now T z xl± as a probability space with a algebra generated by the cylinder 
sets and measure P = ® yeZ d d\x where d/i is a probability measure on T. We introduce 

accordingly the set of i.i.d random variables wJgT labelled by y G Z d and r G I±, with 
common distribution d\i. We construct a diagonal random unitary operator D(a>) on ?^ 
by 

D(a;)|T,y) = e^ r |T,y). (2.5) 

Eventually, we define the one time step random evolution operator 
U U (C) :H^Hby 

U u (C)=B(u>)§, (2.6) 

emphasizing the dependence on the coin matrix C G U (2d) in the notation, which should 
be thought of as a parameter in the problem. This definition, which is the generalization 
of the one considered in [22] to dimensions larger than one, amounts to make the coin 
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matrix that the walker encounters site dependent in a random fashion, so that the 
transition amplitudes are multiplied by a random phase. Explicitely, 

U u (C)\<t,x) = J2^ :+r(T)C r,a\r : x + r(r)). (2.7) 

T6/± 

Note that U U (C) defines an ergodic family of unitary operators: For any a G Z d , let 
T a : T zdx/ ± -> T zdx/ ± defined by (T a cu) y = u y+a , where u y = (uj T y ) T ei ± - Then 

C/T aW (C) = V a U u (C)V- a , (2.8) 

where the unitary operator V a : H — > H is defined by V a |r, y) = |r, y + a), for all r G 7±. 
In particular, VaSV-a = 8, for all a G Z d . 

The proof of our results below requires some regularity of the measure \i and we 
assume the following comfortable hypothesis 

Assumption R: The i.i.d random variables {wJ} 3/e z d ,re/± have an L°° density: 

dfi(6) = l(6)d6, with < I G L°°(T). (2.9) 

We shall work in a "large disorder regime", meaning that the coin matrix C G U(2d) 
of the random walk, viewed as a parameter of the walk, is close to some CV G U(2d), 
such that the corresponding random quantum walk U w (C n ) is strongly localized in the 
following sense: U^(C n ) admits a bloc decomposition with respect to an infinite family of 
deterministic finite dimensional subspaces of H, the direct sum of which coincides with 
H. Consequently, U u (C n ) completely forbids transitions between these subspaces, and 
has pure point spectrum. Moreover, in order to engage in a perturbative analysis, these 
invariant subspaces should allow for the definition of unitary restrictions of U u (C n ), de- 
noted by U^(C n ), to finite volume portions A of the lattice Z d x I±, for which a detailed 
spectral analysis is available. 

As we will show below, matrices from the set of permutation matrices, 

VM = {CV G U(2d) | C v \t) = |tt(t)), Vt G J±, with n G & 2d acting on J ± , } (2.10) 

satisfy the previous requirements. Actually, it will be enough to prove the existence of 
invariant subspaces for 8^ = S(C n <g> I). Note that any permutation 7T G &2d without 
fixed point admits a unique decomposition in cycles with disjoint supports of the form 

vt = 71^2 ...TT n , where 71 j = (Tj, 7t(Tj), n 2 (Tj), 7r m ^' _1 (rj)), (2.11) 

for a set of n < d indices Tj appearing in sequence in the set I± ordered as in (2.1) with 
T\ — 1 and YTj=i m j = 2^- 

As we prove below, our conditions on S w imply conditions on n G &2d'- 
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Lemma 2.1 The operator S> w has pure point spectrum if and only if tt G &2d satisfies 

rrij — l 

£ r^fo)) = 0, Vje{l,...,n}. (2.12) 

s=0 

Remark 2.2 In particular, n has no fixed point. Moreover, rrij is even. 
We can now state our main result: 

Theorem 2.3 Let U U (C) = D(w)S be defined on U = C 2d <g> / 2 (Z d ) by E> = S(C ®T) 
given by (2.4) and 3(uj) by (2.5), under Assumptions R. Let C w G VM. associated with 
n G &2d satisfying (2.12). 

There exists a 5 > such that for all U U (C) with C G U(2d) satisfying \\C — C 7r || C 2d < 
5 we /iawe 

1. Uu{C) has pure point spectrum almost surely; 

2. for any nonnegative p and any ip G H of compact support, it holds for the multi- 
plication operator \X\ p \r, x) = \x\ p \r,x) = ^X^=i max \ x i\j \ T -, X ) almost surely: 

sup|||XTO(C)^ll«<oo; 



3. there exist finite constants g > 0, c > such that for all x, y G Z> d , T,a G I± 



E„ 



sup |(r,x|, f(U u (C)) \<r,y)\ 

/eC(u),|f/||oo<i 



where C(U) is i/ie set o/ complex valued continuous functions on the unit circle U. 

Remark 2.4 Dynamical localization (3) implies exponential localization (2) and spectral 
localization (1), see [17] for details. 

Remark 2.5 Exponential localization (2) says that the state remains localized in a 
neighborhood of the origin, at all times n, with a probability to find the walker in that 
neighborhood arbitrarily close to one. 



Remark 2.6 The results proven in this paper and the method of proof certainly apply 
to more general unitary walks, as those discussed e.g. in [3], when they are randomized 
in such a way that their transition amplitudes are multiplied by i.i.d. random phases. 
In particular, Theorem 2.3 holds true if the elements of the permutation matrix C w are 
supplemented by deterministic phases. 
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Remark 2.7 The assumption that C be close to C„ amounts to saying that the random 
quantum walks we are considering are perturbations of a (trivial) walk in which a quan- 
tum walker moves on finitely many sites only. What makes the result for the perturbed 
case non trivial is the fact that we have a control on the dynamics over infinite times. 

The technical statement which implies Theorem 2.3 reads 

Theorem 2.8 Under the hypotheses of Theorem 2.3, there exist finite constants 5 > 0, 
7 > and c > such that for all C G U{2d) with \\C - C^\\ < 5, all < s < 1/3, all 
x, y G Z d with \x — y\ > 2, all t,<j G I±, all z G C with \z\ G (1/2, 2) \ {1} 7 we have 

E^xKU^-zy 1 \a,y)\ s ) < ce"^. (2.13) 

Indeed, it is well known that the so-called fractional moments estimate (2.13) on the 
matrix elements of the resolvent implies dynamical localization in the self adjoint case, 
as shown by Aizenman and Molchanov in [5]. The adaptation of this statement to a 
fairly general unitary setup is provided in [17], Theorem 3.2, Propositions 3.1 and 3.2. 

Essentially, the first step of the proof of Theorem 2.8 consists in making use of the 
band structure and form of disorder of U W {C) to get an estimate on the squares (or second 
moments) of the matrix elements of the resolvent in terms of the fractional moments 
(2.13), up to an explicit term which diverges as \z\ — > 1. The second step uses a form of 
functional calculus for unitary operators which can be controlled in terms of the second 
moments of the resolvent, see Proposition 5.1 and Paragraph 5.3 of [17], which leads to 
the third statement of Theorem 2.3. The first and second statements are consequences 
of the third one. 

In dimension larger than two, such an estimate has been obtained in some large 
disorder regime or at the band edges of the almost sure spectrum of the ergodic random 
unitary operator under study, see e.g. [19], [17], [7]. When the operator under study is 
off-diagonal, as it is the case for 11^(0), the simpler method of [19] cannot be applied. As 
shown in [7] which addresses the Chalker-Coddington model, a two dimensional random 
unitary model which is off-diagonal, an Aizenman-Molchanov estimate can be derived 
in a large disorder regime from an analysis of the resolvent of finite volume restrictions 
of the operator. While the Chalker-Coddington model differs from the coined walks 
considered here, we can follow the same route. 

The rest of the paper, devoted to the proof of Theorem 2.8 along the lines of [17] 
and [7], is split in several subsections which detail the main steps of the arguments, 
generalizing some aspects along the way. 



3 Proofs 

By convention, we shall repeatedly use the same symbol c to denote constants which 
may vary from line to line but depend only on unessential quantities for our purpose, as 
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will be clear from the context. Also, we shall not distinguish anymore between symbols 
that denote scalar products or norms in different Hilbert spaces. The spectrum of an 
operator A on some Hilbert space is denoted by cr(A). 

3.1 Properties of U^iC^) 

We prove Lemma 2.1, starting with a generalization of Remark 2.2. We show that if § 
is defined by means of an arbitrary C G U(2d) is pure point, then C is off-diagonal. We 
take advantage of the translation invariance of the operator § to make use of Fourier 
transform arguments. 

Proof of Lemma 2.1: For any ip = ^2ip T {y)\ T ,y) £ tl, we define the Fourier 
transform Tip by 

= = Mx)\t) G L 2 (T d ;C 2d ), where 

r£l± 

Mx) = MvK lxy , x e T d c R d . (3.1) 

y&Z d 

Consequently, the Fourier image of S, .FSJ 7-1 on L 2 (T d ; C 2d ), is computed as 

(?&/>)(x) = T Pr'C®\y' + rir')){y'\^ T {y)\r)®\y) (3.2) 

y,y>eZ d ,T,T>eI± 

c T , tT My)e- ix{y+r{T,)) W) 

= e- ixr{T '\C^(x)Ur') = e- ix ^'\C{T^){x)) Tl \r'). 

t'&± t'&± 

Hence, J r §J r_1 acts on L 2 (T d ; C 2d ) as a multiplication operator by the unitary matrix 

C{x) = $(x)C, where $(x) CT)T = 5 atT e~ txr<yT) is diagonal. (3.3) 

Since i 4 C(i) is a trigonometric polynomial, J-'SJ 7-1 is absolutely continuous unless 
the eigenvalues of C(x) are all independent of x G T d . In particular TrC(x) is indepen- 
dent of x, so that d Xj TrC(x) = —ie~ tx ^Cj t j + ie tXj C-j-j = 0, j — 1, 2, • • • , d. Hence 
Q,j — C-j-j = 0, i.e. C is off-diagonal. Specializing to the case C = C w , we get the 
result stated in Remark 2.2. 

We consider now C = C n . Let ttj be a cycle in the decomposition (2.11). Conse- 
quently, for any j G {1, . . . , n}, any Tj G supp ttj 

m,j — l 

SyMHlr,-, £ r^fo))), (3.4) 

s=0 
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which, in Fourier space is equivalent, for any x £ T d , to the eigenvalue equation 

{T^ T-^ix)^) = e-^^^^Tj). (3.5) 



As above, for to be pure point, we need Y^s=o 1 r ( n j( T j)) = 

t 



Conversely, set for any x £ Z d , any Tj £ supp 7ij 

t 

H]: = >i>ai,{y ! ',,)..r • ^V^'fa))), * 0, 1, ... ,m j - l}. (3.6) 



s=0 



Since the sequence of vectors in "Hj corresponds to the successive images of \Tj,x) by 

Syr, "H? is invariant under for any j if Xls^o lr ( 71 " s ( r i)) = 0- Note that T-Q = Wy iff 
\a,y) enl j . Let 

Wj= K\ andW x = 7£. (3.7) 
Tj-esupp je{i,...,n} 

In case 7r consists in one cycle only, i.e. n — 7Ti, "H^ = "H*. The next Lemma, whose 
proof we omit and which follows from the fact that n is a bijection, shows that is 
pure point. H 

Lemma 3.1 Assuming (2.12), the subspaces W x C H, with dim"Hj = m| are such that 
W x JL U\ if j ^ k and § W (W X ) = W x . Moreover, U x with An < dimH x < Ad 2 also 
satisfies § n (H x ) = U x and H = @ xeZd U x - 

Remark 3.2 The last sum is not direct because Ti x fl T-L y ^ in general. 

Remark 3.3 Since D(cu) is diagonal, the random operator U u (C n ) = D(o;)S 7r is reduced 
by all subspaces T-Lx , W x and T-L x as well. 

We shall need the following probabilistic property of the spectrum of the reductions 

U u (C n )\ H -rj. 

Lemma 3.4 There exists Cj < oo such that for any arc A C U of small enough length, 
we have 

P(<r(C/ w (C w )| w? )nA = 0) >l- Cj \A\. (3.8) 
Proof: The eigenvalues of 11^(0^)1^ can be obtained from the property, 

U u (C*) mi \t£ = det(D(w)| w? )][| w? = e^I^, (3.9) 

see (2.11), where 



Hence 

Note that {9l j (uu)} Tj ^ Uj are random variables distributed according to 

6 T J(u) ~ Lj(9)d9 :=(/*/*• ■■l)(6)d6, m r Md convolution. (3.12) 

Consequently, the eigenvalues of [4,(6^)1^ are correlated, and for any arc A C U of 

length |A| < 2tx /rrij we compute the probability that no eigenvalue belongs to A as 
follows: 

mj— 1 

PW^(COki) n A = 0) = 1 - / 9 . , (3-13) 

r=0 Jm i e 3 A 



In particular, Assumption R implies that there exists cj < oo such that Y^r=o f m e ~ 2l7Tr/m i a Lj(0)d9 < 
cAA\, as \A\ ->■ 0. 



Considering the restriction of U U {C^) to = ®j=i,.-. ,nW x , we get from the indepen- 
dence of the 6y(w)'s associated with the invariant orthogonal subspaces which add up 
to H x 

Corollary 3.5 There exists c < oo such that for any arc A CV of small enough length, 
we have 

¥(a(UUC,)\ n J n A = 0) > (1 - c|A|) c . (3.14) 



3.2 Finite Volume Restrictions 

Let us now turn to the definition and analysis of finite volume restrictions of U W {C). In 
order to do that, we consider a generalization of the deterministic quantum walks (2.4) 
which consists in allowing the coin matrix to dependent on the site of the lattice. Let 
C = {C x G U(2d)} xeIi d, and define a unitary operator on H by 

S(C)= PrC y ®\y + r{r)){y\. (3.15) 

yeZ d ,r€l ± 

Remark 3.6 The operator U W (C) can be expressed as §(C(w)) where 

C(u) = {C» G U(2d)} xezd , with (C x (u)) T , a = e^ + " (T) C T , CT , (3.16) 

see (2.7). 

Let A L = {x G Z d | |x| < L} be a box of side length 2L + 1 in Z d , where |x| = 
maxj =1 ... d Given 7r G &2d without fixed point, and L G N, we set 

SL = S < C >- — C « = { C ' f ' ofhe"™ s t' ' + 11 - d < 317 ) 
U%(C) = B(u)S L . (3.18) 
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The finite rank difference U U (C) - U„(C) = is given by 



T^ = 3(co) Pr(C-C n )^\y + r(r))(y\, (3.19) 

\y\e{L,L±l} 

so that, by the Schur condition, see [23] p. 143, for some constant c > uniform in the 
disorder, 

||Tjl| < c\\C-C w \\. (3.20) 

Let 

UAL = H x and U A ° = n x , (3.21) 

{x I \x\<L} {x | \x\>L} 

where the sums are not direct sums. 

Proposition 3.7 The operators E> L and U^{C) leave the subspaces "H Al and T-L a l in- 
variant. 

Proof: It is enough to consider S L , since D(w) is diagonal, and to show the invariance 
of the finite dimensional subspace T-L Al , since S L is unitary. First note that our choice 
(3.17) ensures that @\ x \ = lHx is invariant under S L . Then, because of (3.6), any vector 
\r,y) with \y\ = L belongs to some T-L x , with |x| G {L — 1,L, L + 1} such that l-i x is 
invariant under S L . Hence, vectors \r,y) with \y\ — L — 1 are mapped by S L to vectors 
\r',y') with either \y'\ E {L — 2, L — 1}, that is to vectors of "H Ai or \y'\ = L, which 
belong to some invariant subspace 7i x i, hence \x'\ e {L — 1, L}. Consequently, T-L Al and 
H A % are invariant. H 

We define the restriction of U U {C) to the finite box A L and to its complement A£ as 

U^{C) = Ut{C)\ H , L and f/5(C) = ^(C)|^ A c. (3.22) 

Remark 3.8 The restrictions to boxes centered at arbitrary points v of the lattice, A L + 
v, have exactly the same properties as those stated below, thanks to the ergodicity of 
U U (C), see (2.8). Hence, without loss of generality, we can consider Al only. 

We make use of perturbation theory in C — C n , see (3.20), to analyze the properties 
of the resolvent of the finite volume restriction U Al (C) by comparison with that of 
U Al (C t ). We shall use the following notations for the corresponding resolvents, for 
z U, C E U(2d) 

R*(C, z) = (U* - z)-\ # G {Al, L}. (3.23) 

We need a probabilistic estimate on the position of the eigenvalues of the unitary random 
matrix U AL (C n ). Let dist(z, B) denote the distance between a point z G C and a set 
BcC. ^ 
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Lemma 3.9 For any z U, any LgN, and 77 > snc/i £/ia£ nL d zs small enough, 

F(dist(z,a(U^(C n ))) < < (3.24) 

Proof: This lemma, based on Corollary (3.5), is a generalization to dimension d of 
Proposition 2.2 in [7] which deals with a two dimensional case. We only sketch the 
argument, omitting the details. There are (2L + l) d invariant subspaces % x in "H Al , for 
which (3.14) holds. By independence, for any AcU with \A\ small enough, 

F(a(U£ L (C n )) n A = 0) > (1 - c|A|) c ( 2L+1 ) d . (3.25) 

Now, for z £ U and 77 > small, a ball centered at z of radius 77 intersects U on an arc 
of length of order 77, if not empty. Hence, 

P(dist(z, ( j(C/^(C 7r ))) < 77) < 1 - (1 - c V ) cLd < c V L d , (3.26) 

when 77 and r]L d are small enough. B 

Next, we show that moments of fractional powers of the matrix elements of the 
finite volume resolvent R^ L (C,z) display arbitrary fast polynomial decay in l/L, for 
L-dependent small values of || C — C^W- 

Proposition 3.10 For any s G (0, 1), any a > 0, any p > 1/(1 — s) and any C G U (2d) 
such that \\C -C w \\< i/L 2 ( a P +d ^ +a / s , there exists ac>0 such that 

E^^xlR^iCz^y}^)^^ (3.27) 
for all z^U, all L G N, all |r, x), \a, y) G H 1 ! with \x — y\> 2. 

Remark 3.11 The proof actually yields the more detailed estimate (3.42) under the 
conditions (3.41), which provides two terms: one depending on R^ L {C % ^z) only, via 
Lemma 3. 9, the other containing the effect of the perturbation C — C n . 

Proof: We adapt the strategy of [7] to our d-dimensional model. For a = 0, the estimate 
holds for any C G U(2d), by Theorem 3.1 in [17], a result based on spectral averaging 
which holds for general unitary operators of the form H)(uj)S with D(cu) diagonal as in 
(2.5) and S deterministic, banded and shift invariant. The statement reads in our case: 
for all s G (0, 1) and all C G U(2d), there exists c > such that for all \t,x), \cr,y), and 
all z £ U: 

E(\(r,x\,R#(C,z)\a,y)\ s )<c, (3.28) 

where stands either for the full resolvent R or for R Al , with arbitrary L G N. For 
a > we use the invariance of the subspaces % x by Uu(C n ), first order perturbation 
theory and Lemma 3.9. 
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We introduce the notation 

\t,x) = \a), with a = (ai,a 2 ) = ( T , x ) G Z d x I±. (3.29) 
Let a ~ (3 denote the property a, (3 G Z d x I± are in the same invariant subspace: 

a~fi 4^3U X such that \a),\p)eH x . (3.30) 
Note that \t,x) ~ \&,y) only implies \x — y\ < 1 in Z d , and remark that 

(a\,Rt L (C n ,z)\f3) = 0ifa^p. (3.31) 
Then we use the resolvent identity 

Rt L (C,z) = Rt L (Cn,z) + Rt L {C) {U A -(C n ) - U A ^(C)) R A -(C n ,z), (3.32) 
the fact that only neighboring states are coupled by U U (C), i.e.: 

(a\,U AL (C)\/3) = if \a 2 - &| > 2 (3.33) 

and the estimate 

II U Al (C) - (C.) || <c\\C- ail , (3.34) 

similar to (3.20), to control the the matrix elements of R Al (C,z): For \u) G "H Ai , 
such that \/j,2 — v 2 \ >2,z$.U 

\(^R Al (C,z)\v}\ < (3.35) 
J2 \(^\,R AL (C,z)\a)(a\,(U^C)-U u (C n ))\P)(P\,R AL (C w ,z)\u)\ 

a,@e% d xl± 

dist (2, a [Uui ) ) a ,i3ez d xi ± 

\a 2 -fi 2 \<l 

where c is a numerical constant as the number of sites in the above sum is finite, in- 
dependent of L. To control the denominator, we consider for z G" U and i] > 0, the 
probabilistic event 

G v (z) := G T^ x/± ,dist (z,a(U A L(C n ))) > V ] (3.36) 

and G°(z), its complement. Recall that by Lemma 3.9 we have that if r]L d is small 
enough, P (G^(z)) < ci]L d . Denote by xa the characteristic function of the set A. Now, 
for < s < 1, p > 1/(1 — s) and q such that 1/p + 1/q = 1, we have qs < 1 and by 
Holder's inequality and estimate (3.28) 

(xoswM |</i|,^(C^)H| s ) < (3.37) 

P (G^(z)) * E w ( | <a*|, ifr (C, z) 1 1/) | S9 ) « < c (77L d ) > . 
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By perturbation theory, for some c > 0, ||C — Ctt|| < cr] ensures that 

dist (a (U^(C n )),z) > V ^ dist (a (U^(C)),z) > |. (3.38) 

Thus, by inequality (3.35) we can estimate the complementary part 

XGv{z) (u)\^\,R^(C,z)\u)\ s < (3.39) 



hence, 



s 

c C Ctt 



dist (z,a(u£ L (C n ))) dist (z,a(ut(C))) 



\\C — C n 
< c- 



rf s 



(xcwH K/xl,^^,^^)! 5 ) < c l|g ^f J ' . (3.40) 

Gathering these estimates, we eventually obtain for any < s < 1, p > 1/(1 — s), 
L G N*, and any r] > such that 

small enough and \\C — C w \\ < cr/, (3-41) 

E w ( | (r, x| , ^ (C, z)\a, y) \ s ) < c ^( V L d )^ + ^"f" 11 ^ . (3.42) 

For a > arbitrary, we first fix the scale 7] by i] = Co/L ap+d for some Cq > to deal 
with the first term and then we fix scale of the size of the difference C — C n according 
to HC-CVII = l/L 2 ^ ap+d ^ +a / s . The constant c can be chosen uniform in L e N* so that 
both conditions (3.41) are met, and the proof of the proposition is complete. H 

Finally, we argue that the fractional moment estimates of the finite volume resolvent 
(3.27) imply exponential decay of the fractional moment estimates of the full resolvent 
stated as Theorem 2.8. This quite general property which does not depend on the 
peculiarities of the model is based on geometric resolvent identities which relate the full 
resolvent to its finite volume restriction and on decoupling/resampling arguments which 
yield estimates on expectations of product of fractional powers of resolvents in terms of 
products of expectations of fractional powers of resolvents. The self-adjoint case dealt 
with in [4] was adapted to a unitary setup which is relevant to the present case and to 
that considered in [7], by [17], Section 13. 

We describe the main steps, relying on [17] for the resampling/decoupling estimates, 
slightly improving the argument by taking into account the dependence of the estimates 
in the norm of the difference between U U (C) and U^(C), see (3.19), which is small in 
our case. 
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3.3 Infinite vs. Finite Volume Resolvents 

The geometric resolvent identity relates the resolvents of U W (C), U^ L {C) and Uu L {C). 
We drop the dependences in C G U(2d), oo and z U with l/2<|z|<2to simplify the 
notation below. Definitions (3.19) and (3.22) yield 

U = U L + T L = U Al © U a l + T L , (3.43) 

We keep track of the dependence in t = \\T L \\, where t < c\\C — C n \\, uniformly in L and 
u. We note R = (U — z)^ 1 and 

R L = (U Al © U A °t - z)- 1 = (U Al - zy 1 © (U A i - z)- 1 . (3.44) 

Proposition 3.12 For every s G (0,1/3) there exists a constant c < oo depending on 
s, such that 

E(\{r,0\,R\a,y)\ s ) < ct 2s (l + cfL^ 1 ) (3.45) 

x Yl H\(r,o\,R L \ U )\ s ) mw\,RW,y)\ s ) 

\u)eu \x')en 

L-2<\u 2 \ = L+l L+2<\x' 2 \<L+4 

uniformly in z G" U with 1/2 < \z\ < 2, L G N and y G Z d with \y\ > L + 5. 

Remark 3.13 In the general case where the matrix elements of T L are of order one, 
one can replace the pref actor by cL d ~ x . Since Proposition 3.10 provides arbitrary fast 
polynomial decay in 1/L of the fractional moments of the finite volume resolvent, this 
will turn our to be sufficient for the argument to follow, as in the case [7]. Nevertheless, 
in the sequel, we take advantage of the fact that the matrix elements ofT L are of order 
t < c\\C — CV||, which is small in our case, see (3. 41), in order to optimize on the 
exponents. 

Proof : Recall that T L has non-vanishing matrix-elements only near the boundary 
of the box A L : 

T L \a) = if \a 2 \ < L-2 or \a 2 \ > L + 2 
T L *\a) = if |a 2 | < L-3 or \a 2 \ > L + 3 (3.46) 

and 

(/3|,T L |a) = if \a 2 - /3 2 | > 2. (3.47) 

We do a double-decoupling, once on A L and once on A L+3 . Using the resolvent 
identity twice yields 

R = R L — R L T L R 

= R L - R L T L R L+3 + R L T L RT L+3 R L+?J . (3.48) 
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Since for y G Z d with \y\ > L + 5 and any a,r G I± 

(r,0\,R L \a,y) = (r, 0|, R L T L R L+3 \a, y) = 0, (3.49) 
we get for such vectors, 

(r,0\,R\a,y) = (r, 0|, R L T L RT L+3 R L+3 \a, y), (3.50) 

which is the geometric resolvent identity to be used below. 
Define the boundary of H Al by 

&H L = {(a,0) G (J± x Z d f | (a\,T L \P) ^ 0} (3.51) 
C {(a, fS) | L - 2 < |a 2 | < L + 2, L - 1 < |/3 2 | < L + 1, |a 2 - /3 2 | < 1}. 

Then, an expansion of (3.50) over the boundaries of H Al and "H Ai + 3 gives 

(r,0\,R\a,y)= (3.52) 
^ (r, 0\, R L \u}(u\,T L \u'}(u% R\v}(v\,T L+3 \v')(v'\, R L+3 \a, y). 

(u,u')edH L 
(v,v')£dH L+3 

Taking the power s < 1 and the expectation, we get 

E(\(T,0\,R\a,y}\ s ) < ct 2s (3.53) 

x v(\(T,o\,R L \u)\W\,m\WlR L+3 \*,y)\ s )- 

(v,v')£dH L+3 

Note that (3.44) implies that only the vectors u G H Al and v' G H A< ^+ 3 are to be 
considered in the sum above and give stochastically independent contributions. At this 
point we resort to a resampling argument to decouple the expectations and to the general 
estimate (3.28) to get rid of the full resolvent term. The resampling argument requires 
s G (0, 1/3), which we will assume from now on. The result is stated as Proposition 13.1 
in [17] and it implies: 

For every s G (0, 1/3) there exists a constant c < oo, depending on s, such that 

E(\(T,0\,R\a,y)\ s ) < ct 2s (3.54) 

x £ E(|<T,o|,it^>| s ) H\(vlR L+3 \°,y)\ s ) 

-l Z L+2<\v' 2 \<L+4 

uniformly in z £ U with 1/2 < |2| < 2, L G N and y G Z d with \y\ > L + 5. 

Next we relate R L+3 to R by means of 

R L + 3 = R + R L+3 T L+3 R, (3.55) 
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and expand 

t l+3 = \w)(w\,T L+3 \w')(w'\. (3.56) 

Altogether we get 

mA,R L+ "\°,y)\ s )<n\(vlR\°,y)\ s ) (3.57) 

+ct s Yl mvlR L+3 \w}\ s \(wlR\a,y}\ s ). 

(w,w')edH L +3 

Then, another application of a resampling argument to factorize the expectations 
together with estimate (3.28) eventually yields (3.45), see Proposition 13.2 in [17]. I 

3.4 Iteration 

To obtain the sought for exponential estimate (2.13), we insert the estimate of Proposi- 
tion 3.10 into that of Proposition 3.12 to get 

E(\(r,0\,R\a,y)\ s ) < c{\ + ct s L^ f* ^ - max E(| (x'\, R\a, y)\ s ) 

L a L+2<\x' 2 \<L+4 

= b(L) max E(| (x'\, R\a, y) \ s ), (3.58) 

L+2<\x' 2 \<L+4 

with t < c\\C - C n \\ = \/L^ ap+d ^ +a l s . Taking a > large enough, with < s < 1/3 and 
p > 1/(1 — s), we have 

t s L (d-i) o and b ( L ^ if ^ ^ oo (3.59) 

we can set L e N* large enough so that for L = L , we have b(L ) < 1 and 

E(\(T,0\,R\a,y)\ s ) <b(L ) max E(| (x'\, R\a, y)\ s ). (3.60) 

L +2<\x' 2 \<L +4 

Then it remains to invoke ergodicity in order to control the expectation at the right 
hand side by the same estimates as above, were the distance \x' — y\ can be decreased 
by L. Eventually, we join the points and y in Z d by a sequence of n < c\y\/L boxes 
of side length 2L + 1 along which we iterate n times the estimate (3.60) and use (3.28) 
at the last step to eventually get 

E(\(r,0\,R\a,y}\ s ) < cb(L ) n < ce^ 1 , with 7 = c| ln(6(L ))|/L > 0, (3.61) 

as requested. This finishes the proof of Theorem 2.8, and that of Theorem 2.3. M 

It is of interest to make the size of the perturbation parameter 5 > \\C — C n \\ as 
large as possible in term of the scale L, by optimizing on the set of available parameters 
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a, s,p, given d. Let us explain what we mean by this. The value L > 1 is determined 
by the condition L > L Q implies 

XL) < W < = e < !, where , = {£ ' f ^ - ^> <^ + *»> + 1 ,3,2) 
with 

01 = a(3 + 4ps) + 2(l + d(2s-l)) (3.63) 

2 = a(4 + 6ps) + 3(l + d(2s-l)) 

and 6 > 1. The size of the perturbation is then determined by 

5(L ) = l/ L 2(ap+<i)+a/s = 1/L p = ( e/6o )p//3. (3.54) 

To maximize S(L ), we need to minimize p/(3 over the parameters a,s,p. Without 
performing a complete analysis of the exponent p//3, we can obtain a value, uniform 
in the dimension d, by sending a to infinity and then taking p = (1 +p')/(l — s), 
s — (1 — s')/3, with p', s' > arbitrary small 

p/(3 = 6/5 + o(l). (3.65) 

■ 
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